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and that of two dynamical scalar flelds. It is explicitly demonstrated that the spectral index of 
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I. INTRODUCTION 

The natures on inflation [H-Q in the early universe have been revealed by the recent cosmological observations 
such as the Wilkinson Microwave anisotropy probe (WMAP) ii , the Planck satellite i0 , and the BICEP2 
experiment ii on the quite tiny anisotropy of the cosmic microwave background (CMB) radiation. Owing to the 
release of the recent observational data, in addition to seminal inflation with single scalar field such as new inflation , 
chaotic inflation 0 , natural inflation m , and power-law inflation with the exponential inflaton potential [T^ , novel 
models of single field inflation have been proposed in Refs, [isl [T3|^ (for reviews on more various inflationary models, 
see, e.g., (iTH^ i. 

In addition to inflationary models driven by the scalar filed (i.e., the inflaton field) described above, there have 
been considered the so-called Starobinsky inflation originating from the higher-order curvature term such as 

term^, where R is the Ricci scalar. This model is observationally supported by the Planck results. Such a theory 
can be interpreted as a kind of modified gravity theories including F(R) gravity to account for the late-time cosmic 
acceleration (for reviews on dark energy problem and modified gravity theories, see, for instance, [23l - l3^ i. Various 
inflationary models in modified gravity theories corresponding to extensions of the Starobinsky inflation have been 
explored in Refs. [3^. [35l|. 

In this paper, we investigate inflation in a theory consisting of two scalar fields which non-minimally couple to the 
Ricci scalar and an additional R^ term®. We consider the conformally-invariant two-scalar-field theory in which the 
conformal invariance is broken by adding an R^ term. In particular, we explore the slow-roll inflation in the cases of 
(i) one dynamical scalar field (namely, we set one of two scalar fields a constant) and (ii) two dynamical scalar fields. 
As a consequence, we analyze the spectral index of scalar mode of the density perturbations and the tensor-to-scalar 
ratio and compare the theoretical results with the observational data obtained by the recent Planck satellite and the 
BICEP2 experiment. It is clearly shown that the spectral index and tensor-to-scalar ratio can be compatible with the 
recent Planck results. 

The motivation to propose our theory is to unify inflation in the early universe originating from the R^ term and 
the late-time cosmic acceleration, i.e., the dark energy dominated stage with dark matter. The R^ term is interpreted 
as the contribution from modified gravity, dark energy is described by one of the scalar fields, and dark matter is 


^ Recently, there have also been studied inflationary models with two/multiple scalar fields (or a complex scalar field with two scalar 
degrees of freedom) such as hybrid inflation [T^ and a kind of its extensions Il6l . 

^ Note that the Starobinsky or inflation in the case of no matter is equivalent to non-minimal Higgs inflation considered in Ref. M- 
^ In Ref. [SSl . inflationary cosmology has been studied in a theory with two scalar fields non-minimally coupling to the Ricci scalar. 
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represented by the other scalar field. Furthermore, it seems that multiple field inflation models can fit the Planck 
data better than single field inflation models. Inflationary models with multi-scalar fields [s^ including the so-called 
curvaton scenario [38| have been constructed and the cosmological perturbations in these models have also been 
investigated [s^. We use units of /cb = c = fi = 1 and express the gravitational constant SttGn by 
with the Planck mass of Mpi = = 1.2 x 10^® GeV. 

The organization of the paper is the following. In Sec. II, we explain our model action and derive the gravitational 
field equation and the equations of motions for the scalar fields. We also examine the slow-roll inflation in the case 
of one dynamical scalar field and study the dynamics of the system including the equilibrium points in detail. In 
Sec. Ill, with the conformal transformation |40l |. we explore inflationary cosmology in the Einstein frame. Especially, 
we study inflationary models in the case that the conformal scalar is the dynamical inflaton field and other scalar 
fields are set to be constants. In Sec. IV, we investigate the slow-roll inflation in the case of two dynamical scalar 
fields. Particularly, we consider the resultant spectral index and the tensor-to-scalar ratio in the Jordan frame (i.e., 
the original conformal frame). In Sec. V, we explore the graceful exit from inflation, namely, the instability of the de 
Sitter solution in the present theory. As a demonstration, we concentrate on the case of one dynamical scalar field in 
the Einstein frame, because, as is described in Sec. Ill, in this case the spectral index and the tensor-to-scalar ratio 
can be consistent with the recent Planck results. Finally, conclusions are described in Sec. VI. 


II. TWO SCALAR FIELD MODEL WITH BREAKING THE CONFORMAL INVARIANCE 
A. Model action and its transformation into the canonical form 

Our model action, which consists of two scalar fields </> and u and an E? term, is described as 




{(jP — U^) 
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R+{S/(j)f - {S/uf 




(II.l) 


where g is the determinant of the metric R is the scalar curvature, a{Y 0) is a non-zero constant, s = ±1 is 
a model parameter, V is the covariant derivative, and J{y) is a function of y defined as y = u/(j). The action in 
Eq. (III. II) without the R^ term has been proposed and studied in Refs. (4ll - l48j| . Here and the following, we have taken 
2 k^ = 1. Regarding the property of the action in Eq. (oni), we remark that if there does not exist the R^ term, 
this action is conformally invariant, while that when the R^ term is added, the conformal invariance of this action is 
effectively broken. We also note that our action may be considered to be invariant yet under the restricted conformal 
invariance [d^. By introducing an auxiliary field $, the action in Eq. (III. ID can be rewritten as 


S = 




$ 


+ ^ - ^+2 - {(Y -u^Yj{y) 


The simplest way to transform the action in Eq. (III. 2D into the canonical form is to take the following gauge: 

ci>+±i^^-uY = l. 

With this gauge, the action in Eq. (III. 2D is transformed into 


S= (Yxy/^ I ~ 


2a L 


1 2 


1 - - uY\ + f [(V 0 )^ - (Vu)^] - - uYJiy) \ ■ 


(II.2) 


(II.3) 


(II.4) 


Here, we explicitly explain the properties of the actions described above. The representative feature of the action 
in Eq. (III. ID is that this action is not conformally invariant, while it is scale invariant. To eliminate the non-minimal 
coupling between the scalar fields and the scalar curvature in the action in Eq. (III. 2D . we have used the gauge in 
Eq. (III. 3D . As a result, for the resultant action in Eq. (IH.4D . the scale invariance is broken, although the canonical 
form, i.e., the Einstein-Hilbert term, is recovered. Therefore, it is considered that the form of the gauge in Eq. (III. 2D 
has two aspects. One is to make the action canonical, but the other is to apparently break its scale invariance of the 
action. 

From the action in Eq. (jll.djl . we obtain the gravitational equation: 


1 s 

R^u T ^ (VV i^u) 

1 


1 


2^^“" 1 2a 1 


1 - -u^) - I - (Vu)2] -k - v?fj{y) }> = 0, 


(H.5) 
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where R^v is the Ricci tensor, and the equations of motion for the scalar fields (p a-nd u: 


sap - [i - - uYJ'iy) = o, 


sDm- u 

6a 


^ '-l2 „,2\1 I „,2\ T/„,\ ^ / J,2 


1 - ^i4>^ - - U^)Jiy) - - U^fj'iy) = 0 ) 


( 11 . 6 ) 

(11.7) 


where □ = is the covariant d’Alembertian operator for scalar quantities, and the prime denotes the derivative 

with respect to y. Furthermore, the action in Eq. (jll.411 is represented as 


S = j d^x^g {i? + I [{Vpf - (V«)2] - V{4>, «, J)} . 


Here, V (</>, m, J) is a potential for u and </>, defined as 


y{(p,u,J) = 


1 

^ L 


+iP^-uYJiy)- 


In this case, the gravitational field equation can be expressed by 

1 „ s . 1 

2 

The equations of motion for u and </> also become quite simple as 


R,iy - +1 {yfipyup - y^uv^u) + | } = o ■ 


= 0 , 

sDu — Fu = 0 , 


( 11 . 8 ) 


(11.9) 


(II.IO) 


(11.11) 

( 11 . 12 ) 


with F^ = dVjdp and Vu = dVjdu. 

The flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric ds^ = —dt'^ + a,^{,t)J2i=i 2 3 is assumed, 

where a(t) is the scale factor. The Hubble parameter is given hy H = aj where the dot shows the time derivative. 
In this background space-time, the gravitational field equations read 


2H + 3H^ - -{<p^ - u2) _ iy = 0 . 
4 2 


In addition, the equation of motions for (p and u become 


siP + 3HP) - F^ = 0 , 
s{u + SHii) + Vu = 0 . 


(11.13) 

(11.14) 


(11.15) 

(11.16) 


B. Single dynamical scalar field model 


We consider that p is the inflaton field and the slow-roll inflation occurs. Namely, 


< 


3H(/) 


in the equation of 


motion for p (III.151) and the kinetic energy (1/2) p"^ of p is much smaller than the potential energy F. We also suppose 
that the mass of the other scalar field u is much smaller than the Hubble parameter at the inflationary stage and 
hence the amplitude of u does not vary during inflation. Accordingly, we set u = uo{= constant) at the inflationary 
stage. In this case, it follows from Eq. (111.12^ that 


Vu = -^u 
ba 




- MP^ - u^)J{y) + -{P^- u^fj'iy) = 0 . 


(11.17) 


This has to be satisfied for a value of p{t) and u = uq, although it cannot be true for an arbitrary function J{y). In 
general, we can take mq = 0 and J = J{y^)- The simplest case is mq = 0 and J = 1. For such a case, the scalar field u 
is totally decomposed from equations during inflation. Consequently, the effective potential of the inflaton p is given 

by 

Kff(<^) = ^ (1 - { 2 ^^)^+CP\ (11.18) 
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with C a constant. This expression is appropriate for any function J = J{y^). The number of e-folds during inflation 
is represented as 





24 V^{uA,J) 


(11.19) 


where (()f is the amplitude of (j) at the end of inflation. For the effective potential in Eq. (III.18|) with cf) ^ (j){, we have 




S 

2 


(/)2 ^ 432aC'ln(s2(/,2 + 288aC()>2 
Y s (s2 + 288aC') 


12s) 


3 Inc/) 
s 


( 11 . 20 ) 


Hence, it is seen that there exist the second and third logarithmic correction terms in Eq. (111.201) in comparison with 
the number of e-folds for the quartic inflaton potential as Xcj)'^ with A a constant. 

The slow-roll parameters in the so-called kinematic approach are defined as (for reviews, see, for instance, [Hill) 


e = 


H 


H 

2HH ' 


( 11 . 21 ) 

( 11 . 22 ) 


For the slow-roll inflation, the spectral index rig of the curvature perturbations and the tensor-to-scalar ratio r of the 
density perturbations are described by ©111 

ng = 1 - 6e -b 2?7, (11.23) 

r=16e. (11.24) 

For the slow-roll inflation driven by the effective potential in Eq. (111.181) . the slow-roll parameters are written as 


/ 1 dl4ff((/))A 

" 16</.® 

■ (s2 -b 288aC)(/)2 - 12s ' 

\Ueff((/)) d(l) j 

s 

288aC'(/)4-b (12 - s(/)2)2_ 


2 1 (PVes{(l)) 24 (s2-b 288ctC)(/|2 - 4s 

sEeff(</)) d</)2 “ “T 288aCcj)^ + {12 - 3(1)^)'^ ’ 


(11.25) 

(11.26) 


where we have used Eqs. (111.111) . (111.131) . and (111.141) . With Eq. (111.251) . we find that s has to be negative because 
e > 0. When e < 0, it follows from Eq. (111.211) that H > 0. This means that not the slow-roll inflation but the 
so-called super-inflation can occur. 

If we take the number of e-folds, whose value has to be large such as Ne = 50-60 enough to solve the so-called 
horizon and flatness problem and the values of rig and r suggested by observations, it is apparently regarded that 
there are three equations (111.201) . (111.231) and (111.241) for three independent variables (0, a, C). By solving these 
three equations, we can estimate viable values of our model parameters. However, a and C are incorporated into 
all the equations in the form of a; = aC. Hence, it is necessary to analyze a system of two equations, for example, 
Eqs. (111.231) and (111.201) . whereas r should depend on rig. From this reason, we may try to modify the initial action 
in Eq. (111.11) in the following way. Let us suppose that s is an arbitrary numerical parameter. In this case, we have 
a system of three equations for three variables. We also remark that our potential is well known as a^otential of 
“Spontaneous symmetry breaking inflation (SSBl)” and a viable inflationary model can be constructed [20j, although 
only for positive values of the parameter s. 

According to the Planck 2015 results, rig = 0.968 ± 0.006 (68% CL) 0 and r < 0.11 (95% CL) Q- These 
values are consistent with those obtained by the WMAP satellite [ 1 , The B1CEP2 experiment has suggested r = 
0.20/'/Q g 5 (68% CL) 0, but recently the H-mode polarization of the CMB radiation is considered to the contribution 
from the dust, and not the primordial gravitational waves Q. 

In our model, rig Ri 0.96 can be realized for 50 < < 60 and a wide range of parameter s, but r < 0.11 cannot 

be satished for any values of s. Therefore, we may put s = —1, although it is difficult to produce the value of r 
compatible with the Planck/WMAP data in our model. In Fig. [U we show the values of rig and r in the wide range 
of X and (/> for s = —1. From this figure, we see that we can acquire r k, 0.21 for rig ~ 0.96. In this case, a typical 
value of X is very large. For instance, for N^. = 55 and rig = 0.9603, we hnd </> « 34.7, x k, 2.8 x 10®, and eventually 
obtain r k. 0.212. 

We explore the values of our initial parameters a and C. For clarity, we take Ne = 55. Using the definition of 
X = aC, we find 


Ueff 


1.0 X lOi'^ 


a 


(11.27) 
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Figure 1: Us and r as functions of x and (j> for s = —1. The sheet whose height changes from ~ 0 to ~ 1 denotes ns, while the 
other sheet whose height varies from ~ 0 to ~ 5 shows r. 

This implies that the viable values are C < 3.0 x 10“® and hence a > 1.0 x 10^®. We mention that for such a large 
value of a, the term 1/ (da) in the Friedman equation becomes small, and therefore it could play a role of the effective 
cosmological constant Aeff. This term may lead to the late-time cosmic acceleration. 


C. Equilibrium points 

Next, we investigate equilibrium points in the system. We have the system consisting of two dynamical equations 
(III.Ill) and (III.121) with a constraint equation (III.131) . To examine equilibrium points in this system, we need to rewrite 
the system of two second order differential equations as that of four first order differential equations. It is clear that 
this task is equivalent to explore the shape of the potential, namely, to find its extreme values and study their natures. 
We execute the numerical analysis by using the graphics of the potential for several values of parameters. 

For most of possible shapes of the function J(t/), there is no true minimum point. This is because the value of the 
potential in Eq. (III.91) on the lines = (fP is exactly equal to 1/ (2Qf). This fact is true for such kind of functions as 
J = C, J = C{y‘^ — I)'" with m a constant, J = C cos^(j/^), J = C'exp(?/^), and so on. The typical behavior of a 
potential for such kind of functions is drawn in Fig. [5] Nevertheless, it is possible to construct functions J for which 
the potential in Eq. (III.9|) has a true minimum. For instance, we have 

J{y) = , f ■ (11-28) 

(r -1) 


The typical behaviors of the potential E as a function of (j) and u in Eq. (III.9|) are plotted in Figs. [IHH 














































6 



Figure 2: as a function of 0 and u in Eq. (liLQll near an extreme value for a = 1.0 x 10^®, C = 3.0 x 10 s = —1, and J = C. 


10 - 



Figure 3: 1^ as a function of 4> and u in Eq. (Hm around the minimum for a. = 1.0 x 10^®, C = 3.0 x 10 s = —1, and 
J = C{y^-l)-\ 
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Figure 4: V as a function of (f) and u in Eq. (jll.Qll around a = 1.0 x 10^®, C = 3.0 x 10 s = —1, and J = C{y^ — 1) ^. 


D. Re-collapse and bounce solutions 

We explore the possible re-collapse and bounce solutions^. For such kinds of solutions, the relation H = 0 has to 
be satisfied at some time. At the time, from Eqs. (III.131) and (III.141) . we have 

H=\v. (11.29) 

It is clear that the condition iJ > 0 is necessary for bouncing solutions, whereas the condition H < 0 has to be met 
for re-collapsing solutions. Accordingly, it can be seen from the general form of the potential V in Eq. (III.91) that for 
J > 0, the possibility for bouncing solutions to exist becomes higher because the value of V is defined to be positive. 
On the other hand, if the potential can have negative values, the existence of re-collapsing solutions is possible as 
well. However, in the case that the value of the potential is always positive, it is impossible for re-collapsing solutions 
to exist for, e.g., J given by Eq. (III.29|) with a > 0 and C > 0. 

Related to the bouncing solutions, we mention that the anti-gravity regime in the extended gravity theories with 
the Weyl invariance (4ll - l44 including F{R) gravity [i^ has been examined. 


III. INFLATIONARY COSMOLOGY 

In this section, we reconsider the theory whose action is described by Eq. (jll.lj) and build an inflationary model in 
another way. 


A. Conformal transformation 

The action in Eq. (III. II) is written in the so-called Jordan frame. Instead of taking the gauge in Eq. (III. 31) . we first 
introduce an auxiliary field $ and then make the conformal transformation, i.e., the Weyl re-scaling, of the metric 
from the Jordan frame to the Einstein frame [40j | 

9tJ,v — ^9fiv ) 


^ Recently, cosmological scenarios to avoid the initial singularity in the early universe have proposed in various modified gravity theo¬ 
ries [s^l . 




















where the bar shows the quantities in the Einstein frame. The Ricci scalar is transform as 


R = A-i 


R - SA-^DA + |a-2(VA)2 


By plugging this relation into Eq. (III. 21) and setting 



(III.2) 


(III.3) 


we find® 


. - / \ R - iA-(VA)^ - ^ T ^A(^^ - 


A2(<^2 _ ^2)2 


+ |A [(V0)2 - {Vuf] - A2(<^2 _ J ^ 


(III.4) 


where we have removed the auxiliary field $ with the relation (IIII.3E If we define A as A = e^ with A a scalar field, 
Eq. (IIII.4I) is represented as 


5= l_^e^(^2_^2) 


1 2 


2a 1 12 


+ {(V</-)' - (Vm) 2} - e2"(</>2 - u^fJiy)] . 


(III.5) 


For simplicity, from this point, we will not write the bar over the quantities and operators in the Einstein frame. 
We introduce the following form of the potential 


E(A,(/),u, J) = ^ [l - - u^)] - u^ fj{y ). 

la I 11 J 

The equations of motion for the scalar fields A, 0, and u are given by 


(III.6) 


3nA + - (Vu)2] - Ea = 0 , 

se^n</i + se^V^AV^(/) + = 0, 

se'^Du + se^V^XVuU — Vu = 0. 


(111.7) 

(111.8) 

(111.9) 


Moreover, the gravitational field equation becomes 

R^. - ^Rgf.. + - ^V^AV,A+ |e - |e^ [{Vcj))^ - (Vm)^] + ^(VA)2| = 0. (III.IO) 


In the FLRW background, from Eqs. (IIII.7I) - (IIII.9I) . we acquire 

3A + 9iJA+-u^)+ Ea = 0, (III.ll) 

se^^ + 3se^iJ0 + se^A(^ - E^ = 0 , (III.12) 

se^it + Sse^iJu + se^Au + Ei = 0 . (III.13) 

Furthermore, from Eq. (HIlToIi . we get 

3H^ + -e^((/)2 - u^) - -A^ - Ie = 0, (III.14) 

4 4 2 

2H + 3H^ - - ii^) + -A^ - ^E = 0 . (III.15) 

4 4 2 


^ Clearly, we see that \/—g = 
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B. Inflationary model 


To calculate the observables for inflationary models including the spectral index of curvature perturbations rig and 
the tensor-to-scalar ratio r, it is necessary for two scalar fields to be made constants. If the scalar held A is a constant 
and one of the scalar helds cj) or u plays a role of the inhaton held, we obtain the similar theory to that described in 
the previous section. 

We suppose that the scalar held A is the inhaton held and other two scalar helds are set to be constants as 
(j) = (j)o{= constant) and u = uq{= constant). It should be cautioned that both these scalar helds (j) and u cannot 
be made zero, because in this case, Ves = l/(2a) and hence the spectrum of the curvature perturbations is hat. 
Nevertheless, it is possible for one of these scalar helds (e.g., u) to be taken as zero, while the other has a non-zero 
value. In fact, however, it is impossible to construct any inhationary model with at least (even) one viable parameter 
in this way. This is a reason why the relations (j) = 4’o ^ 0 and u = uo ^ 0 with uq ^ ±^o are forced to be set. 
Therefore, with Eqs. (IIII.12I) and (IIII.13I) . the conditions = 0, 14 = 0 during inhation are expressed as 


Vu 


S 

6a . 
s 

6a - 






1 - - u^) e^uo - - 


ul)uoJ{yo) - 


-ug)V'(yo)^ =0, 

ulfJ'iVo)]- = 0 . 
00 


(111.16) 

(111.17) 


There are the following cases in which these equations are realized: 

• Case 1: The relation K = (s/12)e^(0o ~ '“o) ^ 1 is met. In this case, the hrst terms in the brackets [ ] in 
Eqs. (IIII.16I) and (IIII.17I) may be neglected. Namely, the term with e^ is much smaller than any other terms. 
Accordingly, all the other terms have the multiplier factor e^^ incorporated in the same way, so that this overall 
factor can removed from the other terms. 


• Case 2: The values of the first and second terms in the brackets [ ] may be similar with each other, but the 
overall coefficient term s/ (6a) is sufficiently small. As a consequence, the first terms in Eqs. (IIII.16I) and (IIII.17I) 
are suppressed to be much smaller than all the other terms in these equations. 

Provided that (at least) one of Cases 1 and 2 is realized, namely, the first terms in Eqs. (IIII.16I) and (IIII.17I) are 
negligible, by multiplying Eqs. (IIII.16I) and (IIII.17|) by uq and (po, respectively, and summing them, we acquire® 

J'(2/o) = 0. (III.18) 


Therefore, the function J has to reach its extreme value when the argument becomes y = yo. Moreover, by substituting 
this value into Eqs. (IIII.16I) and (IIII.17I) . we see that in Case I, the relation J{yo) = —s^j (288a) has to be met, whereas 
in Case 2, the relation J(?/o) = 0 has to be realized. Note also that a very wide class of functions may satisfy both 
these conditions. In the following, we study these two cases in more detail. 


1. Case 1 

In Case 1, by combining the value of the function J with the expression of the potential in Eq. (IIII.6I) . we have the 
effective inflaton potential 


I4ff(A) 


I 




(III.19) 


This is some kind of the potential in the modified Higgs inflation model [l^. The number of e-folds in the slow-roll 
regime is expressed as 



(III.20) 


Remind the facts that (j) = (j)Q ^ u = uq ^ and uq 7^ i^Q. 
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with Af the value of A at the end of inflation. If A Af, the representation of Ne{X) for the effective potential in 
Eq. (IIII.191) becomes 


iVe(A) = 


A + 


or by taking into account the assumption that K = (s/12)e^((^Q — Uq) » 1, we may rewrite iVe(A) as 


s(</>o-“o)e^ 

\fA2 „,2\ 


(111.21) 


N,{cj2)^-X. (III.22) 

With Eqs. (III.211) and (III.221) . the slow-roll parameters can be described by using the effective potential as 


1 / 1 dEeff(A)y 

3 \ Ves{X) dX ) ’ 

2 / 1 d2Kff(A)\ 

3 \ Ves{X) dX^ ) ■ 


(111.23) 

(111.24) 


Here, the coefficients of 1/3 in Eq. (IIII.23I) and 2/3 in Eq. (IIII.24I) originate from the non-canonical definition of 
the scalar field A, namely, the coefficient of the kinetic term for A in the action in Eq. (IIII.5I) is 3/2, and not 1/2. 
Eventually, rig and r are described as 


C^e^^ - 60Ce^ -b 108 
3(Ce^ - 6)^ 
IbC^e^'’' 

(Ce-^ — 6)^ ’ 


(111.25) 

(111.26) 


where we have defined (/ as (/ = s(^q —itg). Consequently, we see that rig can be written as a function of r as ng = ng(r). 
However, the values of the parameter lead to those of rig and r, which are compatible with the observations. Indeed, 
for Nf, = 60 and C, = 0.01, we obtain r = 0.004 and ng = 0.9617. In this case, we find K = 0.013 (<C 1). For = 50 
and C = 0.18 (the value of C for W = 50 becomes larger than that for W = 60 by one order of magnitude), we acquire 
r = 0.005 and ng = 0.9568. In this case, we have K = 0.015 (<C 1). As a result, we see that this case is inconsistent, 
because the initial assumption is AT 1, but the consequences suggest AT <C 1. It means that this case cannot be 
realized, in contrast with Case 2, as is shown next. 


2. Case 2 


In Case 2, as already mentioned above, we have the following relations for the function J: J(yo) = 0 and J'{yo) = 0. 
A simple example of such kind of the function J is 


J{y) = c{y -yoY, 9 > 2 , 

where q is a constant. The effective potential is given by 


1 q < .A 




where we have used C = s(</o ~ addition, the expression of W reads 

Ne{4>) = 4^+^- 

Accordingly, we find 

432 - 5C^e2^ - 168Ce^ 


3(Ce^ - 12)2 
64b^e^'’' 


(HI.27) 


(HI.28) 


(HI.29) 


(HI.30) 


3(Ce^ - 12)2 ■ 


r 


(HI.31) 























11 


Thus, it is seen that this theory also has only one parameter, and therefore we obtain rig = ns(r). Through the 
numerical calculations, the value of rig may be set near the value of rig = 0.96 suggested by the observations. In 
fact, in the case that Ng = 60, for C = 10“^, 10“^, 10“^, 10“^, 10“®, and 10“®, we have rig = 0.9652, 0.9641, 
0.9629, 0.9617, 0.9604, and 0.9589, respectively. For all of these cases, the value of r is about r = 0.004. Moreover, 
similarly to the above results, in the case that Ne = 50, for ( = 10“^, 10“^, 10“^, 10“"^, 10“®, and 10“®, we acquire 
rig = 0.9577, 0.9561, 0.9543, 0.9524, 0.9504, and 0.9481, respectively. For all of these cases, the value of r is in the 
range of r = 0.005-0.008. Consequently, in this case, we can obtain rig Ri 0.96 and r < 0.11. These results are 
compatible with the observations of the Planck 2015. 

Furthermore, we examine another aspect of this theory. If Nf^ = 60 with C = 10“^, 10“^, 10“^, 10“^, 10“®, and 
10“®, or if Nf, = 50 with these values of C, namely, the same combinations of values of iVe and C shown above, we 
obtain the estimation as (C/12)e^f « 0.14-0.17. It means that for all the combinations of values of = 50-60 
(during inflation) and those of C, described above, the value of the effective potential I4ff may be estimated as 

0 7 

Feff ~ . (III.32) 

2a 

This relation implies that for sufficiently large values of a, the values of the potential become under the Planck scale. 
This fact is in good agreement with our initial assumption that s/ (6a) <C 1. As a result, it is considered that Case 2 
is quite viable. 


IV. DYNAMICAL TWO SCALAR FIELD MODEL 

In this section, we explore the theory proposed in Sec. II, whose action is given by Eq. (nni, and consider how to 
realize inflation by using both dynamical two scalar fields (f) and vJ. 

The definitions of the slow-roll parameters in Eqs. (III.21|) and (III.221) are used to describe the slow-roll inflation in 
the present theory. In addition, the spectral index rig of the curvature perturbations and the tensor-to-scalar ratio r 
are supposed to be represented as Eqs. (III.241) and (III.231) . respectively. This assumption has been justified in Ref. [^ . 

We consider the FLRW background. We start with the system of the equations (III.Ill) and (III.121) . namely, (III.131) 
and (III.141) in the ELRW space-time. By imposing the standard slow-rolling conditions u <C iiH, cj) and 

77^, the Eriedmann equation and the equations of motion for (f) and u are reduced to be the following 
simple forms: 



9 1 

^ =6^’ 

(IV.l) 


3si7<^ = 14 , 

(IV.2) 


SsHii = —Vu . 

(IV.3) 

First, we derive the number of e-folds by using the initial definition 



^af ntf 



Ne= dlna = Hdt ^ 

Jai Jti 

(IV.4) 

where a\ and Of are the values of the scale factor a(t) at the beginning and end t{ of inflation, respectively. Moreover, 

we also have the following relation 

dV = Vudu + V(i,d(j) = Vuudt + V^j^^dt . 

(IV.5) 

Bv usinff Eos. (IIV.1D-(IIV.3L we sret 

^ 2 - K? ’ 

(IV.6) 


^ It should be noted that in the framework of this theory, infation is realized only for non-interacting scalar fields |54| . whereas in the 
most general case, infation in this theory has not been realized yet, and it is not so clear to analyze the tensor-to-scalar ratio. 
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with (j)f and Uf the amplitudes of (j) a-nd u at the end of inflation, respectively. It follows from the definition in 
Eq. (III.21D that e is described as 


E K? - Vl 

^ 2HV sV^ 

On the other hand, with the definition in Eq. (III.221) and Eq. (IIV.II) . we hnd 


?7 = - 


1 


4IIVV 




+ 2VV 


2 

(^1 - K ?) 


(IV.7) 


(IV. 8 ) 


where in deriving the last equality, we have used Eqs. (IIV.2I) and (IIV.3I) . The expressions (IIV. 6 p - (IIV. 8 l) can be used 
for any choice of the function J, but in the most general case, the calculations would be too cumbersome. In addition, 
the expression (IIV. 6 I 1 leads to only a relation between the scalar fields (/> and u, and not both the values of (j) and 
u simultaneously. Hence, another additional assumption to determine both the values of tp and u (e.g., ~ m) is 

necessary, or one of the values of tp and u may be regarded as a free parameter. 

In the simplest case J = C, it is possible to acquire the values of the main parameters analytically®. In this case, 
from Eq. (IIV. 6 I) . we find 


Ne = ^ - u^) 


432a; 
288x + s 2 


In [(288a; + s^) 


— u^) — 12 s] . 


(IV.9) 


Expressions for e and rj takes the form 


f (288a; + s^){(p‘^ — u^) — 12 s 1 

s ( 288a;((/)^ — + [12 — s ((/)2 — J 

8 ] (288a; + s^) (3(/)^ — 2(p'^u^ + 3u^) — 12 s (^^ — v?) 

_s - m 2 )J 288a; (02 - u'^f + [12 - s (02 _ 


(IV.IO) 

(IV.ll) 


With the expression of in (IIV.9I) . the value of (0^ — u^) can be obtained. It determines the value of e, but that of 
r] cannot be found by using only that of (02 — u'^). Since e(> 0) should be positive, from Eqs. (IIV.9I) and (IIV.IOII . we 
have 02 — ^2 > 0 and s < 0. These relations and x = aC > 0 lead to the positive value of rj. 

We estimate numerical values. For simplicity, we take s = —1. From Eqs. (IIV.9I) and (IIV.IOI) . we get (0^ — u^) and 
X for any interesting values of Ng and r = 16e. Thanks to the structure of the equations, we obtain the value of r 
as r < 0.11. On the other hand, it is seen that rig = where we have defined a new parameter = u^cp"^. 

This relation yields the value of Ug consistent with the observations. Indeed, for = 50 and x = 10®°, we have 
(p'^ — = 2593 and r ~ 0.1. In this case, Ug is represented as 

rig = 3.7 X 10"°C;^ + 0.9815. (IV.12) 

Moreover, for TVe = 60 and x = 10^®, we acquire cp'^ — = 2317 and r ~ 0.11. In this case, rig is expressed as 


rig = 5.0 X 10"°C1^ + 0.9793 . (IV.13) 

If the larger x is, the smaller r becomes, but the minimum value of rig in both Eqs. (IIV.12D and (IIV.13I) increases. 
Equations (IIV.12|1 and (IIV.13P imply that it is impossible to for the pair of rig and r to take their values compatible 
with the observations. The consequence is the same for the different values of s and W- Thus, for the simplest case of 
J = C, it is impossible to realize the Planck analysis with two dynamical scalar fields. However, even in this simplest 
case, the BICEP2 results may be realized. For W = 50 and x = 10°, we obtain <p'^ — = 1267, r ~ 0.202, and 

rig = 3 X 10“®!/^ + 0.9621. Therefore, the case of < 10® is quite viable and consist with the BICEP2 experiment. 

As a result, in this simplest case (of function J), the Planck results cannot be realized. Nevertheless, the model 
with two dynamical scalar fields leads to significantly different results in comparison with the inflationary model with 
the single scalar field. There is another possibility to consider the more viable types of the function J, e.g., Eq. (III.291) . 
In such a case, it is quite difficult to analyze the equations analytically (for example, Ng, is described as an integral 
equation). 


® Even in this simplest case, equations cannot be reduced to the one-field equations because of the kinetic terms. 
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V. GRACEFUL EXIT FROM INFLATION 

In this section, we investigate the graceful exit from inflation, namely, the instability of the de Sitter solution at 
the inflationary stage for the present theory. Especially, we demonstrate the instability of the de Sitter solution in 
the case of one dynamical scalar field in the Einstein frame. In this case, as shown in Sec. Ill B 2, the spectral index 
and the tensor-to-scalar ratio can be compatible with the observations by the Planck satellite. We also examine the 
contribution of an term in the Jordan frame to the instability of the de Sitter solution during inflation. 

We set the perturbations of the Hubble parameter at the inflationary stage as 

iJ = + , |<5(t)|<I (V.l) 

where i/inf(> 0) (= constant) is the Hubble parameter during inflation (whose value is positive), and 5{t) means the 
perturbations from the de Sitter solution. In the present case, only the scalar field A is dynamical, and (j) and u are 
constants. By taking the time derivative of the gravitational field equation (IHI.15() with V = I4ff, where I4ff is given 
by Eq. (IHI.28I) . we have 


2ff + 6HH+hA+:^ 
2 24q; 


l-^eMe^A = 0. 


Moreover, from Eq. (IHI.llI) with V = I4ff (in Eq. (IHI.28I) '). the equation of motion for A reads 

c 


3A + 9HX - 


I2a 


I - —e^ ] e^ = 0. 
12 


(V.2) 


(V.3) 


Since we consider the solution at the inflationary stage in the early universe, we take the limit t —>■ 0. In this limit, 
we could have an approximate solution A ~ In(iJf) for Eq. (IV.31) with the quasi de Sitter solution H Ki 1/ (3t). By 
taking into account this solution, we see that in the limit t —> 0, the third term proportional to A in the left-hand 
side of Eq. (jV.2() is much larger than the fourth term, which is proportional only to A, and therefore the fourth term 
could be neglected. 

To examine the instability of the de Sitter solution, we represent 6{t) as 

S{t) = , (V.4) 


with P a constant. When /3 is positive, the de Sitter solution during inflation becomes unstable, and hence the universe 
can exit from the inflationary stage and then enter the reheating stage. This is because for /3 > 0, the amplitude of 
6{t) increases in time. By substituting the expression of the perturbations in Eq. (IV.II) with Eq. (IV.41) into Eq. (IV.21) 
and using approximate solutions A Ri In {Ht) and H k. 1/ (3t), we find 

Ol 

2iJi„f/32 + f/3 + —Hi, = 0. (V.5) 


The solutions for this equation are obtained as 


P± 


3 (-1 ± yio) TJinf 
2 


(V.6) 


As a consequence, we acquire the positive solution of /3 = /3+ > 0, and thus the universe can exit from the inflationary 
stage. 

We note that even if the other scalar field becomes dynamical, the procedure to examine the instability of the de 
Sitter solution is basically the same as the one demonstrated above. We examine the perturbations of the Hubble 
parameter by using the gravitational field equation with solutions for the equation of motions in terms of two dynamical 
scalar fields. Qualitatively, the form of the solution for the perturbations will be changed, but in principle there can 
exist a solution to represent the property that the de Sitter solution is unstable. 

The contribution of the term in the Jordan frame to the instability of the de Sitter solution is included in the 
scalar field A and its dynamics through the auxiliary field <I> in the Einstein frame. This fact can be seen from the 
action in Eq. (|H.2I) . and the relation (IHI.3I) with A = e^. Accordingly, it is considered that the R? term can be related 
to the graceful exit from inflation, i.e., the instability of the de Sitter solution to describe the slow-roll inflation. 
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VI. CONCLUSIONS 

In the present paper, we have studied inflationary cosmology in a theory where there exist two scalar fields non- 
minimally coupling to the Ricci scalar and an additional term. We have investigated the slow-roll inflation in the 
case of one dynamical scalar field and that of two dynamical scalar fields. We have analyzed the spectral index rig of 
scalar mode of the density perturbations and the tensor-to-scalar ratio r in comparison with the observations of the 
recent Planck and BICEP2 results. 

For the case of single dynamical scalar field in the Jordan frame, if the number of e-folds during inflation is 
50 < iVe < 60, we have ris Ri 0.96 and r = 0(0.1). On the other hand, in the Einstein frame, for the case of one 
dynamical scalar field, we can acquire rig ~ 0.96 and r < 0.11. These are consistent with the Planck 2015 results. 
Furthermore, for the case of two dynamical scalar fields in the Jordan frame, when 50 < W < 60, we obtain rig ~ 0.96 
and r = 0(0.1). As a result, we have found that in the present theory, the spectral index and the tensor-to-scalar 
ratio can be compatible with the recent Planck analysis. We have also shown that in the present theory, the de Sitter 
solution representing the inflationary stage is unstable, and therefore the universe can successfully exit from inflation. 
The term in the Jordan frame is considered to be related to the instability of the de Sitter solution, namely, the 
graceful exit from the inflationary stage. 

As the further developments on the cosmological consistent scenario in the present theory, it is possible to unify 
inflation in the early universe realized by the R^ term and the late-time cosmic acceleration by the dark energy 
component of one of the scalar fields. In this theory, dark matter can also be explained by the other scalar field. To 
construct such a unified scenario is the significant purpose in this work. In addition, the consequences in multiple 
field inflation models seems to be more suitable for the observational data obtained form the Planck satellite than 
single field inflation models. 

It is remarked that the action in Eq. (III. II) may be generalized so that the conformal invariance can be broken by an 
arbitrary function of the F{R) term. Thanks to the additional term in the gravity sector, the novel contributions to 
cosmology yield from the breaking of conformal invariance of the theory. Thus, we have more possibilities to realize 
the unified scenario of inflation, dark energy, and dark matter mentioned above. 

In fact, we may further extend the investigations in this work. We start from the theory with the conformally 
invariance and three scalar fields, and then add an R? term. We explicitly show the expressions of these actions 
below. In Ref. [H^, a theory without an R? term and that with the general potential {(j)^ — v?' — tp^)‘^J{u/(j), u/ip) have 
been explored. The action for the theory with the more specific potential il)'^J{u/(j),u/'>p) and without the R^ term) 
is given by 


5 = 


(0^ — u'^ + e^^) 
6 


i? + (V(^)2 - (Vm) 2 + e(VV')" 




(VLl) 


where e is a constant and ip is the additional third scalar field to two scalar fields (p and u already introduced. If the 
term is added to Eq. (|VI.II) . the representation of the action reads 


S= d^X\/—g 


] ( 


-R^ 


{(j)^ —u'^ + 


6 


R + {Vpf - {Vuf -k e(VV')^ 


u u . 




(VI.2) 


Moreover, when the general form of the potential is included and the signature of e is opposite to that in the action 
in Eq. (IVI.2I) . the action is represented as 


S= 




a 


R^ 


{(jp — V? — ei/'^) 
6 


R+{V(j)f - {S/uf - tiVipf 


u u , 


_(<^^_^^_ei/,^)V(-,-)^ . (VI.3) 


The same analysis as in Secs. II-IV may be adopted to the models in Eqs. (IVI.1I) - (IVI.3I) in which the additional third 
scalar field tp may play a role of the other specie of dark matter, or the additional contribution to the dark energy 
dominated stage or inflation. 
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